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Giant adsorption of microswimmers: duality of shape asymmetry and wall curvature. 
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The effect of shape asymmetry of microswimmers on their adsorption capacity at confining channel 
walls is studied by a simple dumbbell model. For a shape polarity of a forward-swimming cone, 
like the stroke-averaged shape of a sperm, extremely long wall retention times are found, caused 
by a non-vanishing component of the propulsion force pointing steadily into the wall, which grows 
exponentially with the self-propulsion velocity and the shape asymmetry. A direct duality relation 
between shape asymmetry and wall curvature is proposed and verified. Our results are relevant for 
the design microswimmer with controlled wall-adhesion properties. In addition, we confirm that 
pressure in active systems is strongly sensitive to the details of the particle-wall interactions. 
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Boundaries dominate biological process on all scales. 
On the microscopic scale, motile organism may bump 
into various obstactles and boundaries, such as liquid- 
gas or liquid-liquid interfaces, elastic cell membranes or 
solid walls. A universal behavior is the accumulation of 
microswimmers at boundaries [1, 2]. Aside from physico¬ 
chemical effects [3], such as van der Waals forces, two 
main mechanisms have been suggested to explain the 
wall accumulation, hydrodynamic interactions (HI) [4, 5] 
and excluded-volume (or steric) forces [6, 7]. The im¬ 
portance of HI on accumulation is still a subject of de¬ 
bate. However, recent experiment in quasi-2D microflu¬ 
idic channels indicate that surface scattering of sperm 
and Chlamydomonas at lateral boundaries is dominated 
by steric forces with multiple flagellar contacts [8, 9], in 
strong contrast to sperm confined in a 3D channel, where 
HI seems to dominate adhesion [10]. 

Most theoretical studies of simple model swimmers, 
both in bulk and in confinement, have considered so far 
cells with a symmetric body shape, in particular, rods 
[6, 7] or spherical particles [11, 12]. In reality, however, 
cells usually do not exhibit such high symmetry, and 
the stroke-averaged shape of sperm or Chlamydomonas 
rather resembles a forward or a backward swimming cone, 
respectively [8, 9, 13], see also Fig. l(b,c). This raises the 
question how a broken for-aft symmetry of the particle 
shape alters the wall accumulation of cells. 

In order to elucidate shape effects on effective ad¬ 
sorption, we neglect HI and study a generic model of 
an active Brownian dumbbell with unequal bead sizes, 
see Fig. 1(a). Our simulation results show that swim¬ 
mers with a sperm-like (polar) shape exhibit huge wall 
trapping times due to a nonvanishing component of the 
propulsion force directed steadily toward the wall, thus 
resulting in a restricted rotational movement. The trap¬ 
ping times increase exponentially with the shape asym¬ 
metry 9o and the propulsion strength V and could, for 
realistic parameters of Oq and V, exceed trapping times 
due to near-field hydrodynamic forces [5, 14, 15]. In con¬ 


trast, microswimmers with Chlamydomonas-like (antipo- 
lar) shape behave similarly to symmetric rod-like parti¬ 
cles. 

Both in a natural environment and in microfluidic de¬ 
vices [9, 16, 17], microswimmers usually do not swim in 
straight, but rather in curved or branching microchan¬ 
nels. Therefore, the influence of surface curvature on 
accumulation of microswimmers [14, 18-22] is of great 
interest. Based on the analysis of an asymmetric par¬ 
ticle near a flat boundary, we predict a direct duality 
relation between the effect of shape asymmetry and sur¬ 
face curvature on accumulation. For example, a polar 
microswimmer close to a flat wall behaves similarly to 
an apolar particle near a concave surface (e.g. a cavity). 
This is of high relevance for the design microswimmers 
with controlled wall-adhesion properties. 

We model the microswimmer as a self-propelled Brow¬ 
nian dumbbell. The dumbbell consists of two spheres 
with radii Oi and 02 connected by a rigid rod of length 
I, see Fig. 1(a). Its orientation is characterized by a unit 
vector e directed along the axis from sphere 2 to sphere 1. 
The equation of motion for the swimmer’s center (mid¬ 
point between the sphere centers) is then 

r = Ve + a~^Fy, +r], (1) 

where V is the bare propulsion velocity, H is the trans¬ 
lational friction tensor, is the steric force due to 
swimmer-wall interaction and 77 is a random velocity. 
The particle is confined in a channel of height L along 
the z-direction. The sphere a € {1,2} interacts with 
the walls via a screened Goulomb potential Ua [23] with 
a large inverse screening length n and thus resembles 
a hard-sphere. The total dumbbell-wall force is F^, = 
Fi -f F2 with Fq = —VrJJa where Tq, is the position 
of sphere a. The Gaussian white-noise velocity rj obeys 
= 2kBTS^^S{t — t'), where ksT is the ther¬ 
mal energy scale and H = 7|jee-|-7_L(I—ee) is the transla¬ 
tional friction tensor with the friction coefficients 7|| and 
7 j^ for motions parallel and perpendicular to e, respec- 
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FIG. 1. (color online) (a) Sketch of the dnmbbell model of 
asymmetric microswimmers in confinement. The swimmer 
propel along its instantaneons orientation e with velocity V 
in a channel of height L. Left: a antipolar particle (0o > 0); 
right: an polar particle (6^o < 0). (b,c) Experimental images, 
(b) Superimposed phase-contrast micrographs of swimming 
bull sperm; the cell mimics a forward-swimming cone (polar 
shape), (b) A Chlamydomonas alga, confined to quasi-2D 
motion, resembles a backward-swimming triangle (antipolar 
shape). Reprinted with permission from Ref. [13]. 

tively. The orientation evolves according to 

e = (Tu,/7r-b 0 X e, (2) 

where ^ is a Gaussian white-noise vector with 
(l(t)^(t')) = 2Drl5{t — f) and = ksT/^r is the ro¬ 
tational diffusion coefficient. The torque due to the wall 
interaction is Tyj =Ti +T 2 with Ti = (ri — r) x Fi = 
l{e X Fi)/2 and T 2 = —l{e x F2)/2. We solve Eqs. (1) 
and (2) numerically using standard methods [24]. 

Dimensionless numbers characterising the system are 
the Peclet number Pe = V/{lDr), which is the ratio 
of the swimming persistence length V/D^ to the rod 
length I, and the shape asymmetry parameter sin(do) = 
(ai — 02 )//, which is sin(do) < 0 for polar (sperm-like) 
and sin(0o) > 0 for antipolar (Chlamydomonas-lihe) mi¬ 
croswimmers. Wherever appropriate, we choose realistic 
parameters similar to that of Escherichia coli [25]. We 
increase swimming velocity up to Pe = 234; for compari¬ 
son, Escherichia coli achieve Pe « 120 [5], while Chlamy¬ 
domonas and bull sperm reach only Pe « 25 — 50 due to 
the large active rotational diffusion [5, 26]. We consider 
only small asymmetry, |sin(0o)| ~ |^*o| ^ 0.125. Note 
that I sin (0o)l ~ 0-5 for sperm and Chlamydomonas [13], 
see also Fig. l(b,c). Furthermore, we apply a weak con¬ 
finement with channel height L = lOL 

There are universal features in the behavior of an elon¬ 
gated microswimmer confined inside a channel. The 
swimmer performs a persistent random walk within the 
bulk region; when it encounters a wall, a torque, caused 
by steric interactions, leads to approximately parallel 
alignment with the wall; finally, the swimmer can escape 
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FIG. 2. (color online) (a,c) Logarithmically color-coded plots 
of the probability density function (PDF) P[ez, z) at moder¬ 
ate activity Pe = 8. (a) P{ez,z) for a polar {6q < 0) and 
(c) an antipolar swimmer (So > 0). Deterministic trajecto¬ 
ries are indicated by solid lines. A polar swimmer is highly 
localised at {ez,z) = (sin (So), u/2). (b) The corresponding 
orientation-averaged PDF p{z) = P{ez, z) dcz- (d) PDF 
of the orientation P{ez , x « 5) at the threshold of steric in¬ 
teractions (<5 = 1/2 -\- a), see dashed line in (a,c). Swimmers 
leave the wall region increasingly parallel with increasing Pe. 

from the boundary when its orientation, as a result of 
rotational diffusion, points slightly away from the bound¬ 
ary [6, 7]. Generally, swimmers are increasingly localized 
near the wall with increasing activity, as can be seen 
from the density profile p{z) in Fig. 2(b). Shape asym¬ 
metry changes the behavior dramatically. Polar swim¬ 
mers are much more strongly adsorbed as their antipo¬ 
lar counterpart, see Fig. 2(b). Moreover, polar parti¬ 
cles point persistently toward the wall at an angle pre¬ 
scribed by the body shape, see high probability density 
at [cz^z) = (sin (0o), a/2) in Fig. 2(a) and compare to 
the probability density function P{ez,z) of a antipolar 
swimmer in Fig. 2(c). A high probability density near the 
boundary is tantamount to a large wall retention time. 
As can be seen in Fig. 3 and will be discussed in more 
detail below, the retention time of the particles indeed 
increases very rapidly with increasing activity and shape 
polarity. 

In order to understand the giant wall accumulation of 
polar microswimmers, and in particular the huge trap¬ 
ping times, we examine Fqs. (1) and (2) with the focus 
on the two relevant variables, the coordinate along the 
surface normal z and the orientation angle 9, see Fig. 1, 
which implies 

i = F sin (9) -h F^h + ’ 9 -, (3) 

9 =—DrXsm(9) PTw/^t P i- (4) 
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FIG. 3. (color online) (a) Mean time a swimmer remains 
within the wall region (z € [0, 5] or z £ [L—S, L]) as function of 
the shape asymmetry 9o', time is normalized by = l/{2Dr). 

(b) Tu, versus Pe. For polar swimmers Tu, grows exponentially 
with V, while for apolar and antipolar swimmers r™ oc 

(c) Mean time Tb a swimmer remains within the bulk region 
(2 G [5, P — 5]) versus Pe. Tb is independent of 9q and oc 
1/-°'^ for Pe-> L/l = 10. 


The first term on the RHS of Eq. (4) is peculiar for ro¬ 
tational diffusion in 3D [15, 27] and can be neglected 
for Pe ^ 1. The noise obeys {r]{t)r]{t')) = 2DS(t — t') 
with D = ksTh = (D|| -1- 2P)j_)/3 and = 

2Dr5{t—t'). We linearize Eqs. (3) and (4) around the sta¬ 
ble point {z*,9*) = (a/2, do) of a fully absorbed particle, 
and define small perturbation as (dz, 59) = {z—z*,9—9*). 
The equations of motion then reduce to an Ornstein- 
Uhlenbeck process. 


_d 

dt 




(5) 


where we assume small do- Otir aim is to estimate the 
mean escape time Te, i.e., the mean time to reach an 
orientation parallel to the wall, ( 2 , d) = ( 02 , 0 ), from the 
stable position ( 2 , d) = ( 2 *, 9*) by rotational diffusion. In 
order to do so, we reduce the complexity of the problem 
further by neglecting the motion normal to the surface. 


= U=l59^. ( 6 ) 

Here, U is an effective harmonic potential for the orien¬ 
tation angle with spring constant k = —KpjV9o/{d'yr)- 
The escape problem, Eq. (6), is related to a first pas¬ 
sage problem. An exact expression of the mean first- 
passage time from any point along any potential to any 
other point exists [28]; however, due to the complexity 
of this expression, an extraction of the leading contribu¬ 
tions seems unfeasible. A low-noise approximation is the 


Kramers rate theory of crossing a smooth potential bar¬ 
rier [15, 28]; here, we use a heuristic expression for the 
mean escape time over a barrier AU = U{59e) = C/(—do), 
which captures the low- and the high-noise limits [5], 



Note that dp < 0 (dp > 0) for polar (antipolar) particles. 

In the simulations, we measure the mean trapping time 
Tu, as the time during which the swimmer remains within 
the wall region (range of steric interactions: 2 S [0, d] or 
2 G [L — 5, L] with 5 = 1/2 -\- a). In case of polar swim¬ 
mers, this trapping time is an estimate of Tg. As can be 
seen in Fig. 3, we do indeed observe exponential depen¬ 
dencies of Tyj in the low-noise regime, with log (t^,) cx dg 
in Fig. 3(a,b) and log(ru,) oc E in Fig. 3(c). The depen¬ 
dence of Tw on V is nicely consistent with Eq. (7); how¬ 
ever, we observe log(Tu,) oc dg instead of log (re) oc dg, 
as predicted by Eq. (7). Considering the various approx¬ 
imations in the derivation of Eq. (7), like the harmonic 
form of U {59) and the dimension reduction, it is not sur¬ 
prising that we do not obtain a perfect agreement. In 
particular, the parabolic description of U breaks down 
with increasing 59] moreover, 17 is a function of 5z and 
should soften with increasing 5z. 

In case of apolar and antipolar swimmers, which do not 
point persistently toward the wall, the above treatment 
does not apply. Instead, can be deduced from an 
analogy to a semi-flexible polymer adsorbed on a wall 
[6, 7], which predicts a scaling behavior Tw oc see 

Fig. 3(c). This process contributes also to of polar 
particles, but it becomes negligible as compared to the 
escape time over the effective potential barrier AU for 
large polarity and Pe^ 1. 

The mean time ti, a swimmer remains within the bulk 
region (outside the range of steric interactions) is inde¬ 
pendent of dg. As can be seen from Fig. 2(d), particles 
leave the wall region increasingly parallel with increasing 
Pe. In agreement with the semi-flexible polymer anal¬ 
ogy [6], we observe (9) cx at 2 « d in the bal¬ 

listic regime. Thus, with Tb oc L/{Vsin{9)) we obtain 
Tb oc LjV'^l'^ consistently with simulations, see Fig. 3(d). 

With this knowledge, it is easy to interpret a global 
measure of the density distribution, the surface excess (or 


adsorption) T = /g [^( 2 ) — pb\ dz, where pb is the bulk 
density. For a passive hard dumbbell, T = —5/{L — 5) < 
0, while for fully absorbed particles it is T = 1. A rough 
estimate of T is 


T 


fS 

I / p{z) dz — 2pb5 

Jo 


2t„ - Tb25/{L - 25) 


O'T’ 

^ ‘ u, 


■n 


( 8 ) 


Using the scaling of and Tb, see Fig. 3, we obtain 
low- and high-Pe limits, which are consistent with the 
simulation results, as indicated in Fig. 4(a). 

The problem of a microswimmer with shape asymme¬ 
try moving near a planar wall bears a strong resemblance 
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FIG. 4. (color online) (a) Adsorption F versus activity Pe 
for antipolar {9o > 0 ), apolar {9o — 0 ) and polar swimmer 
( 6^0 < 0). Asymptotic estimates of F are indicated by dashed 
lines, r of active Brownian spheres (ABS) is shown for com¬ 
parison. Note that polar swimmers are completely absorbed 
(r = 1) above a critical Pe, by contrast, apolar and antipolar 
swimmers are always partially absorbed, (b) F as function of 
the asymmetry 9o at various Pe. 


with a swimmer moving near a curved wall, see Fig. 5. 
Let us consider first a apolar swimmer in spherical con¬ 
finement of radius R, see Fig. 5(a). In this case, as in 
the case of a polar microswimmer at a planar wall, the 
velocity vector in the stable conformation forms an an¬ 
gle with the tangent plane to the wall at the front bead. 
Thus, in both cases, the microswimmer points toward 
the wall and thus should have very long retention times. 
Secondly, the force of a polar microswimmer towards the 
wall can be partially or fully compensated by a convex 
wall, i.e., for a microswimmer moving at the outer surface 
of a sphere of radius R, see Fig. 5(b). In the case of a full 
compensation, we predict the same accumulation behav¬ 
ior as for an apolar particle at a planar wall. Note that 
an apolar microswimmer would strongly scatter at con¬ 
vex wall. Thus, shape polarity provides the possibility 
for microswimmers to move along curved surfaces! 

Hence, it is obvious to define a generalized asymmetry, 
considering shape asymmetry and wall curvature at once, 
as 00 = 6*0 + ^o’ = (®i ~ o,2)/l + kl/{ 2 R) ior R^ I 
|ai ~ 02!, where k = -1-1 for convex and k = —1 for 
concave boundaries. This allows a unified description 
of asymmetric microswimmers near curved walls, where 
©0 < 0 (00 > 0) implies an exponential grow (algebraic 
decay) of with V. 

We have performed various test in order to verify the 
equivalence of shape asymmetry and surface curvature 
[29]. First, we analyse the behavior of an apolar mi¬ 
croswimmer (00 = 0) close to a surface with a curvature 
ranging from that of a convex to a concave wall. The wall 
retention times r^, as a function of the generalized asym¬ 
metry 00 = Gq+Sq and Pe are shown in Fig. 6(a,b). The 
results closely resemble the corresponding dependencies 
of a asymmetric swimmer near a flat wall in Fig. 3(a,b). 
However, oc V~^ for 0o > 0 in contrast to r^, oc 
for 9q = 0; the oc V~^ behavior reflects a simple 
ballistic escape from convex walls. Further, we have 
tested whether shape polarity can be compensated by 



FIG. 5. (color online) (a) An apolar swimmer confined within 
a spherical cavity of radius R is equivalent to a polar swim¬ 
mer close to a fiat wail provided that the angle 9o ~ —1/{2R), 
between the propulsion force of the apolar swimmer and the 
tangent plane of the cavity at the front bead, equals the asym¬ 
metry 00 ~ (ui — 0 , 2 )/I of the polar particle, (b) A polar par¬ 
ticle near a convex boundary behaves like an apolar swimmer 
close to a flat wall if (ai — 02)/! ~ —1/{2R). 


a negative curvature of the wall. We simulate a polar 
swimmer with fixed 6q = —0.09 close to the surface of a 
sphere at different radii R, i.e., at different wall polarities 
Qwaii _ i/(^2R). The results for as a function of 0o 
and Pe are shown in Fig. 6 (c). Again the similarity with 
the results in Fig. 3(b) is striking. Hence, a symmetric 
rod near a planar wall is equivalent to a polar swimmer 
near a convex boundary provided that shape and wall 
polarity cancel, i.e., 0o « 0 (oi — 02 )/^ « —1/{2R). 

Finally, we want to briefly discuss the the wall pres¬ 
sure in active systems. There are several attempts to 
construct an equation of state for active fluids [30-32]. 
However, this idea has been questioned, because the pres¬ 
sure p, measured as the force exerted on the boundary 
per wall area, should strongly depends on the details of 
the swimmer-wall interaction [33], in contrast to ther¬ 
mal equilibrium. Our results support the latter claim. 
We observe that, in contrast to active Brownian spheres 
where p oc F in the strong confinement limit (Pe ^ L/l) 
[30, 34], p cx for polar swimmers, and p oc for 
symmetric rods and antipolar particles, see Fig. 7. The 
difference between a spherical and a rod-like particle is 
that the latter exert a significant force only during the 
arrival at the wall, and align immediately parallel with 
the boundary due to the steric torque Tyj{ez,z). 

In summary, we have shown that a small shape polarity 
of microswimmers leads to extremly long wall-trapping 
times. The exponential dependence of on V and 0o 
is responsible for a nearly complete adsorption of polar 
particles, in contrast to symmetric rods or antipolar par¬ 
ticles. The pressure p in active system is sensitive to the 
details of the swimmer-wall interaction [33], in particular 
the variation of the asymmetry from antipolar to polar 
change the grow of p with V from sublinear to superlin- 
ear. 
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FIG. 6. (color online) Equivalence of body shape and wall 
curvature. Sq' < 0 corresponds to a motion within a spherical 
cavity and Oq > 0 means swimming close to the outer sur¬ 
face of a sphere. is shown as a function of the generalized 
asymmetry 0o = 0o -|- Oq and Pe. (a,b) An apolar swim¬ 
mer (00 = 0) near curved boundaries {Og ^ 0); the situation 
is illustrated in Fig. 5(a). Filled symbols indicate Tu, of an 
asymmetric swimmer (0o 7^ 0) near a flat wall (8^ = 0). (c) 
A polar swimmer (0o = —0.09) near a convex wall (Q™ > 0); 
the situation is illustrated in Fig. 5(b). In (b) and (c) the 
same power-law decays are indicated by lines. 



FIG. 7. (color online) Active pressure p as a function of Pe 
normalized with the ideal gas pressure At large activity 
pressure grows superlinear for polar swimmers and sublinear 
for symmetric rods and antipolar particles, p of active Brow¬ 
nian spheres (ABS) is shown for comparison [30, 34]. 


We have neglected HI in our analysis to elucidate the 
effects of swimmer shape on wall adsorption. However, 
it is worthwhile to briefly speculate on the interplay of 
both effects. We restrict our discussion to pushers, i.e., 
to micro-organisms and self-propelled particles that gen¬ 
erate thrust behind the body, such as sperm or bacteria. 
The retention time due to HI, corresponds to the 
time to reach an angle 9^^ > 0 starting from parallel, 
where self-propulsion outweighs hydrodynamic attrac¬ 
tion [5]. has been predicted to follow an Arrhenius- 


Kramer-like behavior [5, 15]. For a polar pusher, the 
effective potential for the orientation angle now contains 
two barriers at 0 = 0 (due to shape) and 9 = 9^^ (due 
to HI). Thus, a naive expectation would be that the re¬ 
sulting retention time is a sum of the retention times 
due to shape and HI [35]. However, we expect a syner¬ 
gistic effect, because hydrodynamic interaction will both 
additionally attract the swimmer towards the wall and 
further restrict rotational movement, which increases the 
effective barrier height. In contrast, for antipolar push¬ 
ers, hydrodynamic attraction will favor a configuration 
with orientation pointing away from the wall, and hence 
can be expected to shorten the way to the potential bar¬ 
rier due to HI at 9^^ > 0, thereby ultimately decreasing 
Hl-induced wall adhesion. 

Furthermore, we have shown that wall curvature can 
compensate the effect of shape asymmetry, i.e., asymme¬ 
try and curvature are two sides of the same coin. Thus, 
the combination of both effects can be used in order to 
design artificial microswimmers or microfluidic devices 
for particular tasks. For example, microswimmers could 
be designed which move along surfaces within a porous 
medium [20, 36], while corrugated microfluidic channels 
can be constructed to reduce wall accumulation [9, 22]. 

We thank J. Hu, G. Schiitz, A. Varghese, and R. G. 
Winkler for helpful discussions. We gratefully acknowl¬ 
edge support by the DFG within the priority program 
SPP 1726 “Microswimmers”. 
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